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Abstract 

The w = 1 ^ radiative lifetime of NH (X'^'E~,v = 1, A^ = 0) is determined to be Trad,cxp. = 
37.0±0.5stat I'lo s^ys ™®' corresponding to a transition dipole moment of |/iio| = 0.0540 ^Q^QQ-j^g D. To 
achieve the long observation times necessary for direct time-domain measurement, vibrationally 
excited NH [X'^T,~,v = 1,N = 0) radicals are magnetically trapped using helium buffer-gas 
loading. Simultaneous trapping and lifetime measurement of both the NH(t; = 1,A^ = 0) and 
NH(f = 0, A^ = 0) populations allows for accurate extraction of Trad.exp.- Background helium atoms 
are present during our measurement of Trad,cxp., and the rate constant for helium atom induced 
coUisional quenching of NH(w = 1, A^ = 0) was determined to be fc^=i < 3.9 x 10"^^ cm'^s^^. This 
bound on kv=i yields the quoted systematic uncertainty on Trad^xp.- Using an ab initio dipole 
moment function and an RKR potential, we also determine a theoretical value of 36.99 ms, in 
agreement with our experimental value. Our results provide an independent determination of 
Trad, 10 1 test molccular theory, and furthermore demonstrate the efficacy of buffer-gas loading and 
trapping in determining metastable radiative and collisional lifetimes. 

PACS numbers: 33.70. Ca, 31.15.Ar, 34.50.Ez, 33.80.Ps 



Measurement of molecular abundances in astropliysical, atmospheric, and chemical set- 
tings is critical to the discovery of new physics. A myriad of techniques are used to accom- 
plish this. One in particular, infrared absorption spectroscopy, studies light passing through 
a sample {e.g. an atmosphere or stellar cloud) to determine total absorption. In combina- 
tion with transition strengths measured in the lab, molecular density can be determined. 
However, measurement of transition strengths for rovibrational transitions of radicals is par- 







ticularly difficult [1| . This is generally due to imprecise knowledge of the radical density. One 
approach that removes the need for absolute density calibration is time domain or linewidth 
measurement of the transition lifetime. Despite the intrinsic advantages of this approach, it 
has been seldom used for weak transitions. This is because measurement of lifetimes longer 
than a few milliseconds is difficult; the long measurement time required can be limited by 
coUisional quenching or transit time. A by-product of these experimental limitations has 
)een the increasing importance of ab initio calculations in determining transition strengths 

One "striking" example of the need for improved experimental methods and verification 
of ab initio techniques is the case of imidogen (NH). The NH vibrational radiative lifetime 
Trad is a vital number for astrophysical studies since the Av = 1 transitions of the X^S" 
state fall in the 3 /im atmospheric transmission window. The absolute line strength of the 
0-1 vibration-rotation (V-R) transition (Aio) plays a crucial role in determining the nitrogen 
abundance in cool stars [3|, U, Isl, |6|], which, in turn, is tied to questions of stellar evolution 
and internal mixing. The estimated absolute nitrogen abundance of the sun is based on 
atomic nitrogen (Nl) lines and NH V-R absorption [7|. The lack of precise knowledge of 
AiQ is the dominant source of uncertainty in observed abundances based on V-R lines |8|. 
Despite their importance, measurements of the NH V-R transition strengths are particularly 
lacking [sl. Id, bj . 



Chackerian, Jr. et al. ll| have made the only published laboratory measurements 



of Trad for NH by observing the intensity variation introduced by the vibration-rotation 



interaction (the Herman- Wallis effect 1 2l| ) and combining this with the only measurement of 



the permanent electric dipole moment 



and the published theoretical values 



a 



of NH 



10 



X S ). This experimental determination 
1% [iSj cover a range spanning nearly 50%. 



In this Letter we report a direct time-domain measurement and ab initio theoretical 
calculation of the spontaneous emission lifetime of vibrationally excited NH(f = 1, X = 0) 



radicals. NH {v = 1 and f = 0) is loaded into a magnetic trap using ^He buffer-gas loading. 
The trap lifetime of NH(t' = 1) is measured to be Ttrap,v=i = 33.5 ± 0.5 ms, dominated 
by the radiative lifetime {Trad,io)- Corrections to the measured hfetime due to the finite 
trap evaporation (determined from our measured NH(f = 0) lifetime) and the possibility of 
background helium collisional quenching {ky=i) result in Trad,cxp. = 37.0 ± 0.5stat!^o8'^ys ms, 
in agreement with our calculated value of 36.99 ms. We also measure a tight bound on the 
helium- induced NH(f = 1) collisional quenching rate coefficient of A;„=i < 3.9x 10~^^ cm^s~^. 
In 1998 our group demonstrated the first magnetic trapping of molecules [l6| by buffer- 
gas loading CaH(f = 0). It was discovered that both v = and v = 1 states were coolable 

n 

and CaH(f = 1) persisted in the buffer gas as long as CaH(t> = 0) [ITJ, suggesting that 
the collision-induced quenching of vibrational states in the buffer gas was slow enough to 
utilize buffer-gas loaded molecules to measure vibrational decay lifetimes. In 2005, the 
Meijer group in Berlin used a very different method. Stark-deceleration 18|| and optical 



pumping, to electrostatically trap OH(i; = 1) and perform the first vibrational lifetime 
measurement of cold trapped molecules [19[. Their measurement improved the precision 
of the experimentally determined Einstein A coefficient Aio to better than 4% and was in 
agreement with their ab initio calculated value. Our group recently demonstrated trapping 
of ^E state molecules with the buffer-gas loading of ground-state NH into a magnetic trap 
|20(]. The large splitting between levels of opposite parity in molecules in S states makes 



them unsuitable for Stark deceleration due to nonlinear Stark shifts 18| but well suited 
for buffer-gas loading. We find that unlike CaH, the lifetime of the vibrationally excited 
NH molecules in the trap is significantly shorter than the ground state due to spontaneous 
emission. 

The experimental apparatus and method is essentially the same as that described in Ref. 
|20l |. Briefly, NH radicals are loaded into a magnetic trap from a molecular beam. About 
~ 10% of the NH produced in the beam by dissociating ammonia in a DC glow discharge is 
in the v = 1 state. We do not detect any NH(f = 2). The NH(f = 1) along with NH(f = 0) 
enter the magnetic trap, which is inside a cold buffer-gas cell. The cold cell has a 1 cm 
diameter entrance aperture for loading the NH radicals from the beam so that the buffer 
gas exits through the aperture but is replenished through a fill line. The flow rate of buffer 
gas into the cell is controlled to set the buffer-gas density at a constant level. Low-field 
seeking molecules fall into the 3.9 Tesla magnetic trap and are detected using either laser 
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FIG. 1: Vibrationally-excited trapped molecule density vs. time. The population of the NH 
X'^T,~{v" = 1, N = 0) molecules is monitored using laser-induced fluorescence on the 1-1 band of 
the A ^^ X transition (see inset), and the time profiles are fitted to a single-exponential decay to 
obtain observed lifetimes. The loading process that occurs in the first 50 ms is essentially the same 



as described in 20(], and we begin fitting as least 5 ms after the molecular beam has been turned 



off to ensure rotational and translational themalization. 

induced fluorescence (LIF) or laser absorption on the A *-^ X transition. Laser absorption 
spectroscopy indicates that we typically trap more than 2 x 10* NH(w = 0) molecules. The 
excitation laser is a 672 nm CW dye laser frequency doubled to 336 nm 34|. Thermalization 
of the external and rotational motion of the NH proceeds rapidly, as can be verified from 
their trapped lifetime and the Zeeman-broadened spectrum 20|. 

NH(i; = 0) and NH(f = 1) are detected from the rotational ground state (A^ = 0) on 
the diagonal vibrational lines of the Ri A <-> X transition. Figure [1] shows a measurement 
of the number of trapped NH(f = 1) versus time, i.e. the fluorescence observed due to our 
continuous excitation of the A ^ X transition. Lifetimes are measured by x^-fltting the 
fluorescence to a single exponential and the measured lifetime does not depend on laser power 
in this range. The x^ fits are derived from error bars empirically determined from repeated 
LIF counting measurements and calculated shot noise. Figure |2] shows the measured trap 
lifetime of both NH(i; = 1) and NH(f = 0) vs. buffer-gas density at a temperature of 
615 mK. The vibrational ground state lifetime shows the previously studied dependence on 
buffer-gas density [20|. At low buffer-gas densities, the helium enhances the trap lifetime by 
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FIG. 2: Measured trap lifetime for NH(f = 0) (closed circles, Thoid) and NH(t) = 1) (open circles, 
'^trap,»;=i)- For the NH(t' = 0), the lifetime is set by a combination of the trap depth, the diffusion 
enhancement of the lifetime, and collision-induced Zeeman relaxation. The solid curve is a fit to the 
expected functional form [2Q,]. The lifetime of NH(z; = 1) is limited by the spontaneous emission 
lifetime. 

enforcing diffusive motion of the trapped molecules. As the buffer-gas density is increased 
further, the lifetime is suppressed due to the higher frequency of Zeeman relaxation collisions 
with helium. In contrast, the measured lifetime of the v = \ state (Ttrap,i.=i) shows very little 
buffer-gas density dependence, as is expected for spontaneous radiative decay. 

Figure [3] shows trap lifetimes for NH(f = 0) and NH(f = 1) as the temperature of the 
buffer gas is varied. The NH(w = 0}data show the expected strong dependence of lifetime 
on trap depth, as described in Ref. |20|. In contrast, the NH(f = 1) trap lifetime changes by 
less than 5 ms over this range, which again demonstrates that the lifetime is not dominated 
by either trap depth or collisions. This is again consistent with spontaneous radiative decay. 

The spontaneous radiative lifetime of NH(f = 1) is measured using the following proce- 
dure. We repeatedly measure the lifetime of trapped NH(t' = 1) at 615 mK with buffer-gas 
density riHe = 3.6 x 10^^ cm~^. This measurement is done at the low end of our workable 
buffer-gas density range, resulting in relatively long lifetimes for trapped NH(f = 0), thereby 
reducing the effect of finite trap lifetime on the measured value of Trad,exp.- Data is obtained 
from 1041 trapping cycles and fluorescence intensity time proflles taken over four hours, 
yielding a measured NH(f = 1) trap lifetime of Ttrap,D=i = 33.5 ms with a standard error of 
±0.5 ms. 
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FIG. 3: Measured trap lifetimes for the ground (closed circles) and first vibrationally excited states 
(open circles) of NH vs. buffer-gas temperature. The ti = molecules show the previously studied 



strong dependence on temperature 



20]. The solid curve is a fit to the expected functional form. 



The V = \ state shows very little dependence on buffer-gas temperature. 

In order to extract the 1^0 spontaneous radiative lifetime (rrad,exp.) from the measured 
trap lifetime (Ttrap,i)=i), the finite hold time of the trap (rhoid) as well as additional loss 
mechanisms unique to the vibrationally excited molecules (tiqss) must be taken into account. 
Trad,exp. cau bc Written 



1 



'rad. 



exp. 



-1 _ -1 _ -1 • 

'trap, 11=1 'hold loss 



(1) 



The finite trap hold time effects are measured directly by observing the trap lifetime of 
NH(f = 0). Given that the magnetic moment of NH(f = 1) is identical to NH(f = 0), this 
gives a simple correction of At = +3.5 ±0.1 ms. 

Among the possible mechanisms that would cause population changes in NH(t; = 1), 
background helium collision-induced loss is the only likely mechanism we can identify for 
our system. We will refer to this loss mechanism as vibrational quenching, although it should 
be noted that if the collision-induced Zeeman relaxation for the v = 1 state is different from 
the V = state, it would manifest itself with the same functional dependence on nne as 
direct f = 1 ^ quenching. In order to quantify this possible effect, we fit the data in Fig. 
E] assuming as a free parameter helium collision-induced loss of NH(f = 1) (after correcting 
Ttrap,v=i for the trap hold time). This attempted fit results in a coUisional-quenching rate 



coefficient of 1.1 ± 2.8 x 10~^^ cin^s"^, consistent with zero. This iniphes an upper hmit 
on the colhsional quenching rate coefficient for NH(f = 1) of A;t,=i < 3.9 x 10^^^ cni'^s^^ 
and a consequent systematic correction of the hfetime from this possible effect of between 
—0.8 and +2.0 ms where the negative value would correspond to slower Zeeman relaxation 
for V = 1 than v = 0. This is the dominant systematic uncertainty in our measurement, 



which yields a final value of Trad. 



exp. 



37.0 ±0.5 



f2.0 
Stat _Q gsys 



ms where the second uncertainty 



is the systematic uncertainty (completely due to our measured limit on the collision-induced 
quenching rate coefficieiit). Fractional corrections to the lifetime from blackbody radiation 
2l| . hyperfine splitting 221] and Zeeman shifts 221] are calculated to be negligible (< 10~^). 
We can use this measurement to extract the dipole matrix element for this transition, 
/iio = {4>v=i,N=o\fJ'\4>v=o,N=i) where 0^^7v is the vibrational wavefunction. The spontaneous 
emission lifetime can be written in terms of the vibrational transition dipole moment as 



Trad, 10 



iau^ 



\{(f)v=i,N=o{r)\fi{r)\(j)^=o,N=i{r))[ 



(2) 



where a is the fine-structure constant and uj is the transition frequency. The final rotational 
state for the decay to f = is A^ = 1 due to the parity selection rule. We calculate the 
contribution to the lifetime from the A^ = 3 final state to be negligible (< 10~^). Taking 



UJ 



3092.88 cm ^ from spectroscopic constants in the literature 



m 



2J] gives us a measured 



,+0.0009 



D. 



value of ]/iio] = 0.0540'_gQg-^g 

Comparison of our experimental value of Trad,cxp. to theory is crucial to claim resolution 
of the conflicting published measured and theoretical values. To obtain a state-of-the-art 
calculated value we perform an ab initio calculation of the spontaneous emission lifetime of 
NH(f = 1). The vibrational wave function (j)y=i^N=o{f) was computed by solving the one- 



dimensiona 
method 



Schrodinger equation with the sine-function discrete variable representation 



251 ]. The radial potential was constructed with the RKR method [26|, using the 

=i(r) the centrifugal 



23 



24| . To compute 



Dunham parameters from Ram et al. 
term was added to the potential. The molpro program package [27| was used to compute 
the r-dependent electric dipole moment function ^{r). We used the internally contracted 
multireference single and double excitation configuration interaction (MRCI) method [28| 
and the aug-cc-pV6Z one electron basis. The orbitals were obtained in a complete active 
space self consistent field (CASSCF) calculation. The active space consisted of the la — 
6cr, Itt — 377, and 15 — 25 orbitals. All CASSCF configurations were used as reference 



8 



configurations in the MRCI calculation. The result is a calculated spontaneous emission 
lifetime of 36.99 ms, in excellent agreement with our measured value. 

To get an indication of the accuracy of the ab initio calculation we also computed the 
spectroscopic constants and the lifetime using the MRCI potential instead of the RKR 
potential. This results in a transition frequency that is 0.44% too high, a ground state 
rotational constant that is 0.48% too high and a computed lifetime of 37.29 ms. 

Another indication of the accuracy of the ab initio calculation was obtained by computing 

1 single and 



29 



3G 



3l|. 



the dipole function with the partially spin restricted coupled cluster method wit 
double excitations and a perturbative treatment of the triples (RCCSD(T)) ! 
This is a single reference method and the orbitals were taken from a restricted Hartree 
Fock calculation. Core-excitations were included and the aug-cc-pV6Z basis was used. The 
dipole was obtained in a finite field calculation with electric fields of ±10~^ Eh/aQ. With 
the vibrational wave functions computed for the RKR potential this results in a lifetime of 
37.39 ms. The transition frequency computed with the RCCSD(T) potential is 0.73% too 
high. 

Using the MRCI dipole moment function we also computed the NH(f = 0) permanent 
dipole moment. The result, 1.5246 D, is 10% larger than the 1974 experimental value of 



..38±0.07 D [13[ but it agrees well with the 1996 ab initio value of 1.536 D Paldus and Li 



32l |. who already express doubt about the experimental value. 

While both our measured and calculated values for the transition dipole moment agree, 



they differ significantly from previously published values. Chackerian, Jr. et al. Uj have 



observed the Herman- Wallis effect and combined their spectroscopic data with the 1974 



measured value of the static dipole moment of ground-state NH 13!]. This results in a 



measured value for the transition dipole moment of |/iio| = 0.0648 ± 0.008 D, which they 



present with a calculated ab initio value [ll| of |/iio| = 0.0594 D. These transition dipole 



moments are inconsistent with our result of |/iio| = 0.0540 j'^Q'^Q-^g D. Grevesse et al. provide 
a plot of the calculated dipole moment matrix element (including rotation) vs. A^" that 
gives a value of |/iio| > 0.06 D for A^" = 1 J8|. Cantarella et al. use a series of different 
published potentials and dipole moment functions to calculate dipole matrix elements falling 
in the range of 0.0528 < |/iio| < 0.0618 D [2]. Das et al. calculate and oscillator strength 
corresponding to |/iio| = 0.0526 D [15]. 

Our measured lifetime Trad,exp. = 37.0±0.5stat^Q8sys ms can also be compared to published 



values of the Einstein A coefficient using Aiq = rrad,io- Rosmus and Werner iJ] calculate 



AiQ = 28.7±4 ms, and Dodd et al. lOj calculate A^q = 19.3 ms. Both of these calculations 
neglect rotation, and we find convergence for our calculated value only after inclusion of the 
large basis and active space described above. Despite these discrepancies, our measurement 
is in agreement with our ab initio calculated value and we believe this new value to be the 
most accurate yet reported. 

In conclusion, we have measured the spontaneous emission lifetime of the X^E~(f = 
1, A^ = 0) state of the NH radical in a magnetic trap and performed a new ab initio calcu- 
lation and find these in good agreement. This new value can be used to calibrate molecular 
densities in spectroscopic measurements on the NH 1 — transition. Since previously pub- 
lished values for the transition dipole matrix element are almost all larger than our measured 
value, it is likely that NH densities that have been inferred from V-R absorption on this line 
are lower than the actual density of NH radicals. We estimate that current solar nitrogen 
abundance estimates derived from this transition [7|, |8|, |33|] will require a correction of about 
30%. 

Using buffer-gas loading at low helium densities, we have shown the effect of the buffer gas 
to be small enough to measure the lifetime of rmf^^iofNH] in a magnetic trap with the highest 
precision yet reported. Since the Zeeman relaxation cross section for v = 1 and f = are 
unlikely to be very different, future measurement of the vibrational quenching cross section 
of NH(w = 1) with helium may reduce the systematic error bars of this measurement to 
less than our current statistical uncertainty. Our work demonstrates a new application for 
buffer-gas loading that can benefit the physics, astrophysics, and chemistry communities. 
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helpful discussions, as well as Kevin Strecker for suggesting LD688 laser dye. This work was 
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